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Abstract 

A general and exact critical condition of saddle-node bifurcation is derived in closed form for the 
buck converter. The critical condition is helpful for the converter designers to predict or prevent some 
jump instabilities or coexistence of multiple solutions associated with the saddle-node bifurcation. Some 
previously known critical conditions become special cases in this generalized framework. Given an 
arbitrary control scheme, a systematic procedure is proposed to derive the critical condition for that 
control scheme. 
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I. Introduction 

The buck converter, a type of DC-DC converter, is a nonlinear hybrid system which exhibits both 
continuous and discrete dynamic behavior. It is also a piecewise linear system |[I]> l|2l> IS that the 
dynamics is linear in each time segment separated by the switching actions of the converter. Due to the 
switching actions involved in the operations, the converter is more accurately modeled as a sampled-data 
system IH, lH, |[3l. Instability occurs when there exists a sampled-data pole outside the unit circle in 
the complex plane. 

There are three ways that the sampled-data pole leaves the unit circle, thus causing three typical 
instabilities in DC-DC converters iHl, ||5], ||6]. When the sampled-data pole leaves the unit circle through 
1 in the complex plane, the instability is generally a saddle-node bifurcation (SNB) Q (or pitchfork 
and transcritical bifurcations ||6] which are not typically seen in DC-DC converters). The SNB generally 
involves coexistence of multiple solutions Q, or sudden disappearances or jumps of steady-state solutions 
|[5l . When the pole leaves the unit circle through -1, the instability is a period-doubling bifurcation (PDB) 
|[8l . which generally involves fast-scale subharmonic oscillation. It is a phenomenon where the signal 
ripple has subharmonics. When the pole leaves through a point other than 1 or -1 on the unit circle, the 
instability is a Neimark-Sacker bifurcation (NSB), which generally involves a slow-scale quasi-periodic 
oscillation [31, lH. This paper focuses on SNB. 

Averaged models are generally applied to analyze DC-DC converters. It has been known that most 
averaged models cannot accurately predict the occurrence of PDB (subharmonic oscillation) |[8l. By 
considering the sampling effects and increasing the system dimension, improved averaged models can 
predict the occurrence of PDB ||9], ifTOl . However, even with positive gain or phase margins in the 
improved averaged model, instability may still occur lITTI . Therefore, there is a strong need to obtain the 
exact critical condition of instability. 

Typical control schemes commonly used in DC-DC converters are voltage mode control (VMC) and 
current mode control (CMC) ||T2| . Some SNB critical conditions for particular control schemes have 
been known. For example, in CMC with light loading, an SNB critical condition for the buck converter 
has been reported 1(131 . This paper tries to answer the following questions: 

1) Is there a general closed-form critical condition that directly leads to these known conditions, and 
these known conditions become special cases in the generalized framework? 

2) Some control schemes (such as VMC and CMC) may seem different, but do they share the same 
or similar fonn of critical condition? 

3) Given an arbitrary control scheme, is there a systematic method to derive the critical condition for 
that control scheme? 

The answers to all of these questions will be shown to be affirmative. 

Like the describing function approach ||T41 . harmonic balance ifTSl . |[T6l is a tool to analyze a nonlinear 
system. Based on harmonic balance, an PDB critical condition for the buck converter is obtained in Q, 
lITTl . ifTSl . Following similar approach, in this paper, the exact SNB critical condition is derived. 

Note that, here, the critical condition is expressed in the converter parameter space. For example, 
given a source voltage Vs and load resistance R, an illustrative critical boundary in the parameter space 
(fs, R) is shown in Fig. [T] The critical conditions defines the SNB boundary in the parameter space to 
separate stable and unstable regions. When a converter parameter crosses the critical value, the stability 
(or instability) changes. Critical conditions in closed-forms greatly facilitate the converter design, because 
the quantitative effect of each relevant converter parameter can be clearly seen. 

The remainder of the paper is organized as follows. In Sections |ll] and |llll the operation of the buck 
converter and harmonic balance analysis are presented. A general critical condition in terms of loop gain 
is derived. In Section JVl the critical conditions for typical loop gains are derived. In Sections IVTlVIl 



2 



Unstable region 




Figure 1. An illustrative critical boundary in the parameter space {vs, R). 




Figure 2. A VMC buck converter with a compensator Gc(s). 



the proposed approach is systematically applied to various control schemes. Conclusions are collected in 
Section |VlIl 

II. Brief Review of Buck Converter Operation 

Consider a VMC buck converter shown in Fig. |2] where Vg is the source voltage, Vd is the voltage 
across the diode, Vo is the output voltage, Vr is the reference voltage, y is the compensator output signal, 
and Gc{s) = —y{s)/vo{s) is the compensator transfer function. Denote the ramp amplitude as Vm, and 
the ramp slope as nia = h{d). Denote the cycle period as T, the switching frequency as /s = 1/T and 
let cOs ■= '^nfs. 

Similarly, a CMC buck converter is shown in Fig. |3] where a control signal ic controls the (peak) 
inductor current i^, and y = ic — ih (equivalently, Gc{s) = 1). Note that in Fig. [3j the circuit is arranged 
in a way to have a similar output signal y as in VMC. 

The operation of the converter is as follows. Within a cycle period T, the dynamics is switched between 
two stages. Si and ^2. Switching occurs when the compensator output y intersects with the ramp signal 
h{t) := Vmi^ mod T). One has Vd = Vs in stage 5*1 and f ^ = in stage 5*2. The waveform of Vd{t) is 
a square wave with a duty cycle D. The controlled buck converter is equivalent to a nonlinear switching 
model shown in Fig. |4l Note that Vr adjusts the DC offset and does not affect the loop stability. It can be 
implicitly modeled to adjust the offset of h{t). Let the loop gain be T{s). From Fig.|4l T{s) = VsG{s)/Vm, 
and the model in Fig. |4] can be normaUzed as shown in Fig. |5l 
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Figure 3. A CMC buck converter. 
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Figure 4. An equivalent nonlinear switching model. 



III. General Critical Condition Derived From Harmonic Balance Analysis 

The signal path in the control loop has two parts: from y to through a nonlinear PWM modulator, 
and from Vd to y through a linear transfer function G{s) := Gp{s)Gc{s), where Gp{s) is the power stage 
transfer function and its representation depends on what signal is fed to the compensator 

In VMC, Vo is fed to the compensator. Let p = R/{R + Rc). For R^ = 0, p = I. From 111] p. 470], 
the power stage Vd-to-Vo transfer function is 

In CMC, ii is fed to the compensator (with Gc{s) = 1). From |[T2l p. 470], the power stage Vd-to-ii 
transfer function is 

CZs _)_ 1 

^^ + {^ + RcC)s + l 

Let x^{t) be the T-periodic solution of the converter. Let y^{t) be the corresponding T-periodic 
compensator output signal. The intersection of h{t) with y^{t) determines the duty cycle and hence the 
waveform of V(i{t). By "balancing" the equation y^{t) = h{t) (written in Fourier series form) at the 
switching instants, conditions for existence of periodic solutions and SNB can be derived. 

Let d = DT. In steady state as shown in Fig. [6l Vd{t) is T-periodic and it can be represented by 
Fourier series (harmonics) 

^—jnujsd 

Vd{t) = Vs Cne^""^"^ where c„ = (3) 
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Figure 5. An equivalent normalized nonlinear switching model. 




Figure 6. Illustrative signals of h{t) and Vd{t). 



From Figs. Fig. |2] and |4l one has 

oo 

y''{t) = G,{0)vr-Vs Cne^'^^^'GijnUs) (4) 

n=—oo 

Since intersects with h{t) at t = d = DT, one has y^{d) = h{d). Using (O, one has 

oo 

Ge(O)^;,- ^(e^""='-l)G(jna;,) = Mt^) (5) 

This is an equation of d. When SNB occurs, only a single solution exists, the curve of y^{d) is tangent 
to h{d). The SNB is also called a tangent bifurcation |[6l. One has y^{d) = h{d) = ma, which leads to 

oo 

Yl e^'"''''G{jnUs) = ma (6) 

n=— oo 

Note that, using © and the fact that y^{d) = h{d), 

Gc(0)f^ - h{d) 



(V) 



which shows Vs as a function of d = DT. 



Let Re denote taking the real part of a complex number. Note that uogd = 2ttD, e^^^=''^G{jnuJs) + 
e-j^^^d.G{-jnuJs) = 2Re[eJ"'^='^G(jna;,)], and T{s) = G{s)vs/Vn, = G{s)vs/maT. Then, by simple 
algebra, Q leads to the following theorem. 
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Theorem: Consider a closed-loop buck converter shown in Fig. |4] or Fig. \5\ The critical condition of 
SNB is 

- ^ReE e^-^-^GOn^.)] - ^ = (8) 

n=l 

For Vm 7^ 0, ([S]) is expressed in terms of loop gain as 

oo 

- 2Re[ J] eJ'2'^"^r(jnL^,)] = r(0) + 1 (9) 

n=l 

It will be shown later that these critical conditions can be expressed in a closed-form related to the 
hyperbolic function csch. Note that ^ is an expression of convenience. For Vm = 0, the loop gain 
VsG{s) /Vm would be infinite. In that case, the equivalent critical condition ([8]l is used. Since both ([8]l 
and ^ are exact critical conditions, they can be used as benchmarks to determine the accuracy of other 
critical conditions. 

Generally the power stage has an order of at least two (associated with L and C). If the compensator 
has an order of three, the loop gain T(s) has an order of five. However, the compensator is generally 
designed to cancel some poles or zeros of the power stage, and the order of T{s) is reduced. Note 
that the right side of dill involves T(0), but the left side of ^ involves only frequency higher than ajg. 
The shape of the Bode plot for the frequency smaller than is irrelevant. One can use a simplified 
high-frequency form of T{s) for stability analysis. By decomposing the loop gain into partial fractions, 
many new closed-form critical conditions in terms of the converter parameters can be obtained. 

L-plot and S-plot. Define an "F-transform" of T(s) as 

oo 

T[T{s)] := -2Re[ J] e^-2"-^r(jna;,)] (10) 

n=l 

Note that T\T{s)] is a function of many converter parameters. Here, this function is called an L-plot, 
denoted as £. For example, let it be a function of D and it becomes C{D). Then, the critical condition 
(|9l) becomes C := T\T{s)\ = T(0) + 1. The critical condition itself does not tell which side of the critical 
boundary will be the stable region. Generally for a converter, the region with £ < r(0) + 1 is stable. 
Also, define an S-plot as 5 = ma£/(T(0) + 1), which is also the left side of Q or ([8]l. Then, the critical 
condition ([H) becomes S = ma, and the region with S < ma is stable. The S-plot is a useful design tool 
because it shows the minimum ramp slope to stabilize the converter. The F-transforms of typical loop 
gain functions are presented next. 

IV. "F-Transforms" of Typical Loop Gain Functions 

The loop gain is generally designed to have sufficient gain and phase margins. For a stable converter, 
the phase of T{ju]) is less than 180° at the crossover frequency uj^- One can focus only those loop gains 
of first or second orders. 

The loop gain can be further decomposed into a combination of partial fractions. Only partial fractions 
of first orders are considered. Similar analysis can be applied to partial fractions of second orders. 

Let LOp and coz be the pole and zero of T{s). Let p = LOp/cOs and z = lOz/cOs- The F-transform of the 
fraction 1/(1 + oOp) will be the building block to derive the F-transforms of other loop gains. 
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Table I 

F-TRANSFORM OF TYPICAL LOOP GAIN T{s). 



Case T{s) T[T{s)] (note: p — ujp/uja and z = ujz/lOs) 



Ci 


1 


— a(D,p) = — (ao(_D) - ai(D)p + c(D,p)) 


Ca 


1 

s 




Ca 


1 

l + s/a;p 


pa(D,p) 


C4 


l + s/a;p 


p(l-fHAp) 


Cs 


1 


^(ai(Z?)p - c{D,p))= ^(aoiD) - a{D,p)) 


Ce 


1 




Ct 


l + s/u:^ 




Cg 




^(fao(75) - (f - - ciD,p))) 


Cg 


s^(l + s/up) 


^(Ra,{D) + il~^Jc{D,p)) 



From (fTOl l. 

°o ^j2kTTD' 

y- — 

= l/p-7re^P(^-2^)csch(7r^j) (11) 
:= a{D,p) 

where the proof of ([TTI i can be obtained by looking up from a handbook of mathematical formulas 
or checked by a simple computer program, and csch is a hyperbolic function. Using Taylor series 

oo 

expansion, let a{D,p) = ^ {—l)'^ak{D)p'^ . Using the L'Hospital's rule, one has ao(^) = vr(2D — 1) 

fc=0 

and ai{D) = tt'^{2D'^ - 2D + 1/3). A plot of a{D,p) is shown in Fig. Ul The right straight-line edge 
in Fig. His a{D,0) = ao{D) = -k{2D - 1). 

oo 

Let the correction term be c{D,p) = J2 {—^Y OLk{D)p^ ■ One has 

k=2 

c{D,p) = a{D,p)-ao{D) + ai{D)p (12) 

A plot of c{D,p) is shown in Fig. [8] From Fig. [H the correction term c{D,p) is significant only if 
p > 0.2 (equivalently, Up > 0.2oOs)- For p < 0.2, c{D,p) can be ignored. 

The F-transforms of other loop gains are shown in Table U For each case Ci-Cg, the F-transform can 
be derived or proved (by simple algebra) in three ways. First, follow the definition of the F-transform as 
in ( fTOl ). Second, by decomposing T{s) into a combination of fractions 1/s, l/(s + iOp), etc, T[T{s)] is 
a combination of T[l/{s + tOp)], etc. For example, 

• Use of Ci and C2 leads to C5 and Cg; 

• Use of C2 and Cq leads to C7; 

• Use of Ci, C2 and Cq leads to Cg; and 

. Use of Ci and the fact that = leads to C4. 



(jJsTl 1 = -2Re 

S + U}p 



Figure 7. Plot of a{D,p) for case Ci. 




Figure 8. Plot of the correction term c{D,p). 



Third, each case is a special/general case of other cases. For example, by setting ojp — )■ 0, Ci leads to C2; 
C5 leads to Cq; and Cg leads to C7. By setting uo^ 00, C7 leads to Cq. One sees that Ci is a building 
block for other cases because they have similar terms as Ci. 
Remarks: 

(a) All of the transforms in Table U are exact. No approximation is assumed. 

(b) There is no correction term c{D,p) for C2, Cq, or C7. All other cases have a correction term c{D,p), 
which is small and can be ignored if p < 0.2 as discussed above. 

The systematic procedure to derive the critical conditions for two control schemes is discussed next. 
For other control schemes, the procedure is similar. Without loss of generality, assume ESR Rc = 0, then 
p=l. 



8 




Figure 9. A buck converter under multi-loop state feedback. 



V. Multi-Loop State Feedback: Cases C2 and 

Consider a buck converter under multi-loop state feedback, y = Vr — kii^ — k^Vo, as shown in Fig. |9] 
From O and ©, one has T(0) = ki/R + k^ and 



which is a combination of cases C2 and Cq. 

Let K := 2L/RT which is a dimensionless parameter related to the loading as defined in |[T2l . From 
(|9l) and Table IH the critical condition is 



Tr^~^ + rn 2 > = Tr^~5 + fc^j + 1 (14) 



which can be arranged as 



^iD-—)^^(-l^^^ ^)=». (15) 

where the left side of ([TST i is also an approximate S-plot. For ^ 12LC, (fTSl) becomes 

v^kj , „ K +1^ v^k„ 

—)-^ = ma (16) 

which can be rearranged as 

„ K + 1 Lrria Lk,, 

" = — + ^^fir <"> 

Example. Consider a buck converter under multi-loop state feedback ||5l p. 232]. In lH, a digital control is 
used, but it can be approximated as a multi-loop analog control, with y{t) = Vr — kiii — kyVo intersecting 
with h{t) to determine the switching actions. The converter parameters are T = 400 fis, L = 20 mH, 
C = A7 fiF, R = 22 n, ki = 2.1435, k^ = -0.1383, Vr = 0.2152 V, and Vm = 1- 

Using Vs as the bifurcation parameter, the bifurcation diagram is shown in Fig. [TOl and one sees that 
SNB occurs at Vg = 20 and D = 0.7. Generally in SNB, there is a hysteretic loop as shown in Fig. (TO] 
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Figure 10. Bifurcation diagram showing stable (solid) and unstable (dashed) solutions. SNB occurs at Vs = 20, D = 0.7, and 

Vo = VsD = 14. 



3000 
2500 

Q. 
CO 

2000 
1500 




Duty cycle 

Figure 11. The S-plot (solid line) intersects with nia = 2500 (dashed line) at D = 0.7 where SNB occurs. 



In the figure, the upper soUd Une is for the operation when the switch is always on (hence D = 1), and 
the dashed line and the lower solid line are for unstable and stable T-periodic solutions respectively with 
a duty cycle less than 1. For Vs between 19.25 V and 20 V, there are three solutions: a stable T-periodic 
solution, an unstable T-periodic solution, and the third (stable) solution for the switch being always on 
(hence D = 1). When the converter operates with the stable T-periodic solution and Vg is increased a 
little above 20 V, the output voltage will jump up from 14 V to 20 V. Similarly, when the converter 
operates with D = 1 and Vs is decreased a little below 19.25 V, the output voltage will jump down from 
19.25 V to 10 V. The jumping up and down forms a hysteretic loop. 

Using (|7]l, the S-plot based on (fTSl ) (which is almost identical to the exact S-plot based on ([8]l) is 
shown in Fig. [TT] which indicates that SNB occurs exactly ai D = 0.7. The S-plot also shows that if 
the ramp slope is higher than 2898, the converter is stabilized. Also, from ([TtT i and one also has 
D = 0.7 when SNB occurs. □ 
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VI. CMC: CASEC2 

In CMC shown in Fig. ^ y = ic — ii- CMC can be considered as a special case of multi-loop control 
with ky = and ki = 1. One has T{s) « Vg/VmLs, which is of case €2- From ([TTl l with k^ = and 
ki = 1, the critical condition is 

K + 1 Lrua 

D = ^- + (18) 

It agrees with |[T3l which derives the same condition based on either sampled-data analysis or steady-state 
analysis. It is interesting to note that three (sampled-data, steady-state and harmonic balance) completely 
different analyses lead to the same condition. This also corroborates the accuracy of the derived critical 
condition ( fTSl ). 

Without the compensation ramp {nia = 0), the critical point is D = {K + l)/2. 



VII. Conclusion 

A general and exact SNB critical condition ^ in terms of the loop gain is derived. The critical 
condition is applicable to a general nonlinear switching system represented in Fig. |5] Therefore, it is also 
applicable to other similar nonlinear control systems. The effects of different parameters (such as Vs, R, 
and the ramp slope nia) on the instability can be clearly seen. The derived critical conditions agreed with 
the past research results by sampled-data analysis or by time-domain simulations. The critical conditions 
also show the required ramp slope to stabilize the converter. The critical conditions have many different 
forms and lead to different plots. In the L-plot, the critical condition is C = T[T{s)] = 1. In the S-plot, 
the critical condition is S = ma- The derived critical conditions are helpful for the converter designers 
to predict or prevent some jump instabilities or coexistence of multiple solutions associated with the 
saddle-node bifurcation. 

The questions asked in the Introduction are answered: 

1) The previously known critical conditions, such as ([TST l. become special cases in this generalized 
framework. 

2) A typical critical condition is a weighted combination of three terms: ao{D) = 7r{2D — l), ai{D) = 
7r^(2D^ — 2D + 1/3), and a correction term c{D,p). Different control schemes have similar form 
of critical condition if the loop gain function is of the same case in Table Jl 

3) Given an arbitrary control scheme, a systematic procedure is proposed to derive the critical condition 
for that control scheme. First, approximate the loop gain T{s) (for frequency higher than w^). Then, 
from Table U one can readily obtain the critical condition in terms of converter parameters. 

Given the closed-form condition in terms of converter parameters, one knows the quantitative effect of 
each parameter on the instability. One can make the S-plot (such as Fig. [TT]) as a function of a parameter 
of interest, its intersection with the ramp slope nia determines the stable operating range of that parameter. 
The S-plot also shows the required ramp slope to stabilize the converter. 

This paper focuses on the converter operated at a fixed switching frequency. Similar analysis can be 
applied to the converter with variable frequency control, such as constant on-time control. This paper 
also focuses on the SNB. Similar analysis can be also applied for PDB (subharmonic oscillation). The 
F-transforms in Table I still apply, but with a different a{D,p) and a minor modification. The results are 
reported separately. 
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